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Abstract 

We compute the finite size spectrum for the spin 1/2 XXZ chain with twisted 
boundary conditions, for anisotropy in the regime < 7 < 7r/2, and arbitrary twist 
9. The string hypothesis is employed for treating complex excitations. The Bethe 
Ansatz equtions are solved within a coupled non-linear integral equation approach, 
with one equation for each type of string. The root-of- unity quantum group invariant 
periodic chain reduces to the XXZx/ 2 chain with a set of twist boundary conditions 
(71-/7 £ ZZ, 9 an integer multiple of 7). For this model, the restricted Hilbert space 
corresponds to an unitary conformal field theory, and we recover all primary states 
in the Kac table in terms of states with specific twist and strings. 
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1 Introduction 



Apart from their relevance to real lattice systems with restricted dimensionality, 
spin chains are interesting from a field theoretical point of view as integrable lattice 
regularizations of 1+1 dimensional field theories 0. 

The continuum limit of the field theory is the thermodynamic limit of the spin 
chain, where the number of spins in the chain goes to infinity If the thermodynamic 
limit is accompanied by taking the lattice spacing to zero, keeping the length of the 
chain constant, we get a description of the continuum field theory in finite volume. 
At criticality, one expects conformal invariance, with the conformal weights read out 
via the finite size corrections of the energy and the momentum [Q, [5[ Q. 

In the class of theories with trigonometric i?-matrices, massive theories (sine- 
Gordon/massive Thirring models) are usually associated to spin-1/2 anisotropic 
Heissenberg chains (XXZ1/2 chains) with inhomogenieties 0, whereas massless 
theories are related to homogeneous chains. The conformal properties of the latter 
are well known, the central charge c = 1, and conformal weights are ones of the 
gaussian model @, 0, §]• 

Central charges c < 1 can be reached by considering chains with some special 
fixed, or twisted boundary conditions, related to critical Potts and Ashkin- Teller 
chains || H [I0|, |IT|. The correct interpretation of this lowering of c in the case of open 



chains was found by Pasquier and Saleur [12]. The key feature is the invariance of the 



spin chain dynamics under the universal quantum enveloping algebra U q (sl2)- (The 
quantum deformation parameter is given by the spin chain anisotropy as q = e 17 ). 
For open chains, the boundary conditions can be chosen to ensure U q {sl2) invariance. 
This, however, is not sufficient to alter the central charge, which is a bulk effect. The 
lowering of c comes about for a root-of-unity q (rational 7/vr), when the quantum 
group invariance allows for a self-consistent truncation of the Hilbert space. By 
restricting the Hilbert space to type-II (a.k.a "good") representations of the quantum 



group, the central charge is lowered from c = 1 to those of minimal models ||12j| . 
Furthermore, analyzing low lying excited states shows that the conformal weights 
of primary operators form a one-parameter subset of the Kac table. 

For closed chains, quantum group invariance is a more delicate matter |13, 14]. 
A topological interaction of the Wess-Zumino-Witten type has to be introduced [13] , 
which leads to non-local terms in the Hamiltonian. On the level of Bethe Ansatz 
equations, the quantum group invariant closed chain can be viewed as a collection 
of XXZ^ chains with a set of twisted boundary conditions commensurate with 
the anisotropy. Again, quantum group reduction leads to a lowering of the central 
charge to c < 1 f|TJ], |T5[. The collection of twist boundary conditions is the exact 
analogue to the #-vacua of SOS and sine-Gordon models, in which the quantum 
group restriction leads to RSOS or RSG models, respectively. The quantum group 
invariant chain can be considered as the critical/ultraviolet limits of these models 
(see e.g. Ref. JT§). 

The main motivation for this work is the interpretation found in Reference [17], 
of the conformal properties of twisted XXZ1/2 chains, in terms of a lagrangean 
conformally invariant field theory with c < 1, namely the Liouville model with 
imaginary coupling. If the imaginary Liouville model on a circle is discretized 
preserving the integrable structure, it is equivalent to the quantum group invariant 



periodic XXZ1/2 chain on the level of Bethe Ansatz; the Liouville model can be 
mapped to a set of XXZ1/2 chains with twisted boundary conditions. For root of 
unity couplings, a quantum group reduction can be performed on the Hilbert space. 
The central charges of minimal models are reproduced, confirming the dependence 
of central charge on Liouville coupling, known from canonical quantization. 

Furthermore, based on Karowski's treatment of Potts models |H| , primary states 
furnishing the full Kac tables of unitary minimal models were conjectured to arise 
from a specific set of Bethe Ansatz states. More exactly, primary states should arise 
from BA states with one string excitation above the sea of real BA rapidities, and 
no extra holes, i.e. vanishing total spin S of the chain. 

This leads us to reconsider the finite size analysis of XXZ1/2 chains with twisted 
boundary conditions. In the extensive literature on the subject |[9], |7|, p], |9|, |I0 



20, BTL E2L [2^, |4], |25|, 15], the main stress has been on investigations of pure 



hole excitations. Very little has been done to analyze complex excitations. String 
excitations with a restricted set of couplings and twists were investigated in fl8|j , 
and complex pair excitations were considered in [Q. Apart from these, most work 
on complex excitations has been numerical ||, [24], f25| . 

We shall use a non-linear integral equation (NLIE) method [ 261, 27] to analyze 
the finite size spectrum. This method is more powerful than the Euler-MacLaurin 
methods mostly used in the literature, and it allows for extracting the scaling 
behaviour from the Bethe Ansatz equations in a rather straight forward manner. We 
follow the Destri-de Vega (DdV) approach [^7[ , extending it to complex excitations 
with a twist. In Reference |Z5|| , the DdV equations for the sine-Gordon model were 
generalized to complex excitations described as wide and close pairs. Here, we make 
the crucial departure from the methods of f28| in that the string picture of the 
excitations is used. This is because, as pointed out above, the primary states are 
most transparently described in this language. 

The plan of the paper will be the following: In Section 2, we formulate the 
albegraic Bethe Ansatz for the XXZx/ 2 chain with twisted boundary conditions. 
This we do for generic anisotropy in the antiferro magnetic regime < 7 < n/2 
and generic twist angle 9. We describe the complex excitations as strings. The 
DdV equations are derived for the case where we have positive parity strings of one 
specified length k in addition to a bulk of real Bethe Ansatz roots. To facilitate 
the analysis, we restrict the range of string lengths to obey k'y < tt. In Section 3, 
we compute the finite size corrections to the energies and momenta of these states. 
We specialize to states where no holes punctuate the distribution of real roots. In 
Section 4, we specialize to the Hilbert space pertinent for lattice Liouville model, 
i.e. the restricted quantum group invariant chain. We show how the quantum group 
reduction lowers the central charge. Analyzing excited states with only one longer 
string in the sector S z = 0, we show that these states are primary states of the 
theory, furnishing the whole Kac table. 

The main point of this paper will be the fact that taking Bethe Ansatz 
configurations with one positive parity string, we are able to retrieve the spectrum 
of unitary minimal conformal field theories using a NLIE method. The two integers 
in the ensuing Kac table are the string length k and the twist k. Aside from the 
technical advantage gained from this way of doing things, the hope is that this 
method is easier to generalize to the less explored case of c > 25 of relevance to 
Liouville theory with real coupling. 



2 Twisted Bethe Ansatz equations for strings 



We consider a spin 1/2 XXZ chain with 2N sites, with anisotropy 7, and boundary 
conditions twisted by the angle 6: 



+ _ -fci0 _± . 2 _ 2 

°2AT+1 _ e °1 1 °2iV+l — °1 



The corresponding hamiltonian reads: 

2N 

Hxxz = 5 £ a t a j+i + Vjrf+i + 2 cos ^ a j a j+i ( 2 ) 

where the antiferromagnetic XXX limit corresponds to 7 = and the free Fermion 
point is at 7 = | . 

We formulate the theory in the standard language of quantum inverse scattering 
|[29|| . The L- matrix, satisfying the fundamental commutation relation RLL = LLR, 
is 

. / sinh(A + i7<7*) a~ sinh i7 



-^n,a(A) 



sinh i7 



(3) 



\ cr+ sinh i7 sinh(A — i7cr*) 

where the operators a^, are defined by the commutation relations of the quantum 
group U q (sl2). The L-matrix is normalized so that for A = i 7/2 it degenerates to 
a permutation operator between the auxiliary and quantum spaces. Using this L- 
matrix, we can write the transfer matrix of the chain, and build the eigenstates of 
the transfer matrix using the algebraic Bethe Ansatz. 

We will focus our attention on those eigenstates of the transfer matrix that can 
be described in terms of real roots Xj, j = 1, ...ni, as well as strings of complex 
roots Xj\ j = 1, . . . ,ni. A (positive parity) /-string A^ is the collection of / roots 
of the form A + irry , n = —(I — l)/2, . . . , (I — l)/2, with A real. Accordingly, in 
compact notations real roots are frequently designated as one-strings. 

We shall only consider positive parity strings. Negative parity strings can be 
treated within the same context but they are irrelevant for our purposes since they 
do not lead to primary conformal states in the thermodynamic limit. 

The number of /-strings is rii, and the total number of roots is J2i I ni = N — S, 
where S is the number of Bethe roots lacking from the maximal number of roots 
N, i.e. the ^-component of the total spin of the chain. In order to simplify the 
underlying analysis, we consider only the antiferromagnetic regime 7 < tt/2, which 
covers all unitary theories. We also adopt a technical restriction on the range of the 
string length, so that < it. This restriction, although not required in the most 
general case of anisotropy, shows up rather naturally when we consider the roots 
of unity case of the XXZ model pertaining to unitary conformal field theories, as a 
consequence of the underlying quantum group structure. 

(k) 

The Bethe Ansatz equations solved by the rapidities A) describing a Bethe state 

are 



e -2ikd 



sinh(Af } + ik-f/2 
sinh (x\ k) - i&7/2 



2N 

= II IlMAf-Af), V(M), (4) 

l, n,^0 j=l 



where the string-string scattering matrix is 



g — (fe+ pj 2 1 sinh(A + irrvy) sinh (A + i(m + 1)7) 
H m =\k i\/2 sinh ~ im ^) sinn - i (771+1)7) 

To solve these equations, it is beneficial to write them in a logarithmic form. For 
this, we define the phase function 

, / \ \ , sinh (im7 + A) 

m A = ilog , ^>™>0 6 

sinh (im7 — A) 

0o(A) = 0. 

For real argument A, the phase function m is a continuous monotonic function. We 
choose the branch of the logarithm so that it is antisymmetric with respect to the 
origin. Using the phase function, we further define the momentum function 

Pj(A) = I/2 (A) , (7) 

the string-string scattering phase 

(fc+O/2-1 

*w(A)= £ [0 m (A)+0 m+ i(A)] , (8) 

m=\k-l\/2 



and the counting functions for /-strings 



zi(X) = 2N Pl (X) - J2 *ik(X ~ Af ) + 219 . (9) 

By defining counting functions for strings, all the analysis is reduced to dealing 
with functions of real variables, a considerable simplification over the usual method 
of treating complex NLIE's. In this respect, the string hypothesis is necessary, 
(although probably not inescapable), and marks a departure from the standard 
Destri-de Vega treatement. Of course, we must keep in mind that this is only 
justified when we are at zero temperature. 

In terms of the counting functions, the Bethe equations (|j) become 

Z^X?) = 27Tif , (10) 

where the quantum numbers I® encode the choice of branch of the logarithm. They 
are integers or half integers when ni is even or odd, respectively. These quantum 
numbers should all be distinct for a specified string length, otherwise the Bethe 
Ansatz equations are easily seen to lead to a infinite repulsion of the quasiparticles 
corresponding to equal quantum numbers. This is what is usually referred to as the 
fermionic character of Bethe states. 

The basic assumption one makes in treating the Bethe equations for real 
rapidities A is that the monotonicity of 2Np± is suflicent to make the counting 
function Z\ overall monotonic. This allows one to find the spectrum of the quantum 
numbers Ij . Indeed, assuming monotonicity, zi/2tc takes all the (half) integer 
values between z%{— oo)/27r and ^i(oo)/27r exactly once. The only freedom we are left 
with in determining the spectrum of the I:s lies in the choice of the overall branch of 
the counting function z±, which is irrelevant. Any possible non-monotonicity would 



induce additional holes close to the ends of the rapidity distributions |p8[ . These 
don't effect the quantities calculated here. As opposed to z\, the higher string 
counting functions z k are generically non-monotonous. 

For a configuration characterized by the higher string occupation numbers n k , 
the asymptotic values of the counting functions are 



Zi(±oo) = =brnj ± 2S(n-lj) ± 2vr^(fc - l)n k + 219 . 

k>l 



(ID 



If n k 7^ for some k, it is clear that there are more vacancies between Z\{— oo) 
and Zi(oo) than those occupied by the Bethe Ansatz roots. These unoccupied 
vacancies are called holes, and they correspond to zeros of the eigenvalue of the 
transfer matrix (as opposed to Bethe Ansatz roots, which correspond to poles with 
vanishing residue). The number of holes in the distribution of real roots is 



h 



k>l 



l)n k + 2S 
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(12) 



Here |_°J denotes the integer part of o. The integer part terms vanish for small 
values of the total spin S. The salient feature of Relation ([12]) is that a /c-string 
gives rise to 2{k — 1) holes. 

The antiferromagnetic ground state is given by the configuration where there 
are no complex roots, and the number of real roots is maximal, i.e. n\ = N; n k = 
0, k > 1; S — 0. From Equation (|12]), it is easily seen that there are no holes; in 
the ground state all vacancies are filled. 

In the thermodynamical limit, low- lying exitations above the ground state are 
characterized by configurations where some higher strings are exited, as well as 
possibly some extra holes (S might differ from zero), but the only macroscopic 
occupation number is the one of real roots, i.e. only rii ~ N. 



Now we make use of the Destri - de Vega method |27[ to derive an integral 
equation equivalent to the Bethe equations. The derivative of the function 



A(X) = log{l + (-1 



(13) 



acts as a density of real roots and holes in contour integrals; it has first order poles 
at all allowed 2ir x (half) integer values of zi(A), with residue one. Accordingly, for 
a function / analytic within the contour, 



I /(A) A' = £ f(Xj) + J2 /(A* 

' ■ 1 m:holes 



(14) 



j : roots 



The integration path T surrounds the roots and holes counterclockwise. Stretching 
T to ±oo ± ie, one gets the coupled equations for the counting functions of real 
roots and /c-strings 

h n k 

z l (X)=2Na(X) + 2\m{G ll * e A)(X)+ £ F n (A-/i m ) - Y / F lk (X-Q + C x (15) 



m=l 



z k (X)=2lm{G lk * e A)(X) + ^ F lk (X - fx m ) - ^F^X-Q) + C k . (16) 

m=l j=l 



These equations generalize the results of Ref. [27] to encompass strings and holes. [] 
We have denoted the rapidities of the holes with \i m , m — 1, . . . , h, and the rapidities 
of the /c-strings with Q, j = 1, . . . n^. The symbol * e denotes a convolution along a 
infinitesimally shifted real line, 



(G * £ A) (A) = / ^G(A -fi-ie) A(fi + ie) . 



(17) 



The inhomogeneity function a, arising from the action of an inverse convolution 
on the momentum function, is 

a(A) = (1 + $n) _1 *pi = arctansinh^A . (18) 

Similarly, the functions F, giving the hole and string contributions, are defined as 

F u = (l + ^iJ- 1 *^ . (19) 

The measures of the convolutions have 2tt denumerators, as in flTT|), and 1(A) = 
2tt5(X). The kernels G of the integral equations are given by differentiation, 
Gu(X) = J--Fii(A), and the effect of the fc-string on itself is encoded by the function 

Fkk = §kk - Gik * $ifc. 

It is remarkable that the equation doesn't have any inhomogeneity arising 
from the momentum functions. This follows from the identity pk = G\k * pi, which 
can be proved by Fourier transforming. 

Finally, the constant terms in the equations are 



Ck 



7T 9 7T(7T — 27) 



7T 



7 2(tt-7) 
(2k - 2)tt 9 7r(7r - kj) 
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7i 
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7T 
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7T 
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7T 



7T 
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7 



7T 



(20) 



(21) 



These have the contributions of of the twist angle 9, renormalized by the integration 
over the inverse convolution The integer part terms arise from boundary 

terms in partial integrations. 

For computing the finite size corrections, one takes iV large but not infinite. From 
the functional form of the source term, 2Na, one can see that the Fermi points, i.e. 
the limits of the distribution of real rapidities, are close to A = ±logiV. To treat 
the excitations close to the Fermi points, one defines the so called kink counting 
functions, 



rXA) = 
kW = 



lim (zi(A ± I log4iV) t Nit) 
Urn ^(A±^log4iV) . 



(22) 
(23) 



The configurations of holes and strings that can be treated exactly in this limit are 
such that all holes and strings are either close to one of the Fermi points, or well 
away from both. Here, we shall assume that every hole and string is in the vicinity 
of one of the Fermi points. 



lr The generalization to pure hole excitations was performed in Ref. jli 



We parametrize the hole rapidities as 

/i m = /x±±^log4iV, m = l,...,h ± . (24) 

The number of holes in the vicinity of the left (right) Fermi point is h + (h~). 
The total number of holes is h = h + + h~. Similarly, the fc-string rapidities are 
Q = ± - log4iV, and are the numbers of strings close to each Fermi point. 
At one Fermi point, the contributions of the holes and strings situated close to the 
opposite Fermi point are just given by F(±oo). 

The only non-trivial scaling behaviour comes from the function a. Indeed, 
subtracting the bulk contribution, one has for values of A close to the Fermi points 

\im (2Na(\ ± £ log4iV) t A%) = =f e T ^ A . (25) 

We get the following coupled non-linear integral equation for the kink counting 
functions: 

z^(X) = ftj(A) + 21m(G ll * e A±)(\) , l = l,k, (26) 

where we have used the scaled density integral A± = \og(l + (— l)™ 1 e lz± and the 
inhomogeneity functions 

m) = TSi,i e"7 A + £F u (A-/i£) - + C±- (27) 

m=l j=l 

The constant terms carry the renormalized twist 9, as well as the information from 
the opposite Fermi point: 

Cf = d + h^F u (±oo) - n*F lk (±oc) (28) 



For extracting energy and momentum eigenvalues from the scaled NLIEs 
one needs to know the asymptotic behaviour of the scaled functions. At ±oo we 
have 

^(±oo) = ±7in k (k-2) ±tvSt 2S-1 + 29 . (29) 
At =Foo, the kink equation (|26| ) implies 

zt (Too) = TzfiToo) = too . (30) 

Due to the infinitesimal e in ImA(A + ie) = Im e A, the value of this expression 
always between — 7r/2 and n/2. In the limit e — > 0, it is a sawtooth function, with 
a step-function jump of — it at each root or hole position. A small finite e makes 
the function analytic. The asymptotic values of the imaginary parts of the scaled 
density integrals are thus 



Im e A ± (±oo) = Im e A(±oo) = 9 =F S7 + n 



\- e -±sl 

2 ix ix 



9 ± T S 1 , (31) 



It should be noticed that close to ±00 the imaginary part operation can glide through 
the Gu to act only on A, up to a correction that vanishes like eloge. That is, the 
integrals in the convolutions of Equation fl2"6|) are taken along the real axis, so that 
the kernels Gu are convoluted with the regularized sawtooth function Im e A 

In the opposite asymptotic regime, deep in the bulk of filled real rapidity states, 
the small but finite e lets ^'(^oo) dominate in the expansion of A±(too + ie), and 



Im e A ± (=Foo) = 



(32) 



3 The Finite Size Spectrum 



At the external rapidity A close to i7/2, the dominating contribution to the 
eigenvalue of the transfer matrix is given by 



A(A; {A?>}) 



sinh(A + i7/2) 
sinh i7 



1 2N n, • LnC) 



n n 



sinh(A t; - A + i(/ + 1)7/2) 



n^ t=i sinh(Af - A - i (/ - 1)7/2) ' 

(33) 

We expect local integrals of motion at the point A = i7/2, where the L-matrix (£J) 
degenerates to a permutation operator. The corresponding momentum and energy 
eigenvalues of a Bethe state are 



P{{Xf}) = ilogA(i T /2;{Af}) = (TrSm + e + SwW) mod 2 (34) 

1,3 

^({Af}) = i^logA(A;{Af}) = 2tfl cot 7 - I E^°) ■ (35) 



\ ; - 7T 7T , • 



Using the logarithmic Bethe equations flTOD , we express the momentum in terms of 
the quantum numbers /: 

p ({^}) = (*En, + ^E J f ) mod 2 ■ (36) 

'J 

To compute the energy eigenvalues, we once more employ the contour integral 
method (13]) for the sum over rapidities of real roots A«- . Doing this, we get 

£({/U;{C;}) = 2Ar^Lt 7 -/^^(Ay(A)^ 

+ li me f 00 — cr'(A)A'(A) + ^E^Vm) (37) 

The dependence of the string rapidities Q is only indirect, coming through the Bethe 
equations in the NLIE form, when evaluating the last two terms. 

The 1/N corrections to the energy come from the parts of the second integral 
in ( |3"TD that are close to the Fermi points A = ±2 log4A r , as well as from the third 
term. From Equation (|TBp we have cr'(X) = - (2 cosh ^A) _1 , which is exponentially 
peaked around the origin. From the asymptotic analysis of Im e A, we however know 
that for large N, lm e A' vanishes as a double exponential when one goes in to the 
bulk from the Fermi points ±^ log 4 A. Thus the main contributions to the second 
integral in (|37|) for large N come from the vicinity of these points. In this regime, a' 



takes an exponential form, (r'(/i±^ log4A) « e^T^ . Also, the hole rapidities in 
the last term of ( p7|) are taken to be close to one of the Fermi points, as in Equation 

After these considerations, the finite size corrections to the energy are given by 

E FSC = ^ (£+ + £-) , (38) 
where the contributions of the right and left ends of the rapidity spectrum are 



£ ± =lmJ^ e T >A' ± + V ^™ ■ (39) 
J ix , 

m=l 



Using the scaled NLIEs (p6|) , these can further be reduced to 



m=l 3=1 J 71 

(40) 

where the Q:s are defined in Equation (|27D . 

The terms in (|40| ) involving the constants C*i and of Equation ( |28| ) can be 
evaluated using the asymptotic calues of the F-functions: Fn(±oo) = ±n =F \^zr 

and Fi fc (±oo) = ±A;tt t (k - 1)^. 

The values of the counting functions at holes and roots are given by 2-7T times the 
(half) integer quantum numbers I . To evaluate the sums over the hole quantum 
numbers in ([40]), we take a distribution where all real roots are in the bulk, so that 
no holes intermeddle with the roots (Fig 1.). Thus the quantum numbers 1^ for the 
positive rapidity holes get the h + largest (half)integer values of ^J t "(A)/27r, and the 
quantum numbers I~ for the negative holes get the h~ smallest (half) integer values 
of zr(A)/27r: 



J£ = ±l((k-2)n k + S-l) 



9 -±sl 

7T 7T 



=F 1 ± m, m = 1, . . . . (41) 



The quantum numbers I k,± of the strings may take any (half) integer values allowed 
by the asymptotics of the counting functions (pTTj) up to non-monotonicities that are 
consistent with the non-linear integral equations (Effi). 

Finally, the integral in the last term in ( |40|) can be evaluated by the standard 
dilogarithm trick, see e.g. [27, ^8|. For the densities A± and the inhomogenities Q^ 1 
satisfying equations of the form fl26|) , and the boundary conditions (|31| , |32|) for A±, 
we have 

y±,f.A a ( . .\ , 1 ^ ,n _ o.A2 

-OO 



Im. 



dfj, Q ± '(fi)A±([i) 



± 



2 7i — 7 



7T 



12 



(42) 



Combining the contributions of the various terms in Eq. (flUD, the finite size 
energy of the state described above evaluates to 



E 



FSC 



2N { 7T - 7 
+ 27r( kAhAn k + Sn k (k - 1 



6> + nAh - ii(k - l)An k 



+ (tt- 7 )S 2 



TT 

6 



;i + 2 A;-2A; 2 ) + ^/f )+ -^/ 



r(*)- 



j'=i 



j'=i 



where the antisymmetry of the fc-string distribution is An* 
antisymmetry of the hole distribution is 



n. 



n 



k ■ 
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In ix 



2 7T 7T 
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(44) 
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bulk of N—S—kn^ real roots 



Figure 1: The distribution of real roots and holes in a state with one k -string. 



With the distribution of real roots as described, it is an easy task to compute 
the momentum (0). This turns out to be 



P = n(N - S -kn k + l- Ah) mod2vr + — (kn k Ah + s(-+Ah) + £ifM 

J ~ l (45) 

The twisted ground state is the state with no strings nor holes. For zero external 
magnetic field, it lies in the sector 5 = 0, and has vanishing momentum. The ground 
state energy is 



2N \7T-7 

The central charge is c = 1, as the well-known formula of Refs. |2], [3] connecting the 
ground state energy with the central charge gets corrected in the presence of non- 
trivial boundary conditions. The central charge is a bulk property, which cannot be 
changed by boundary effects. 

For the subsets of parameter space investigated earlier, the computed excited 
state energies and momenta agree with results in the literature. For the pure hole 



states (no complex excitations), they coincide with the ones in Ref. |j21|| . Regarding 



states with complex excitations, the two-string result agrees with Ref. [p2| . For 



generic strings, our results agree with the ones of Ref. fll8| , for the discrete set of 
couplings (7 = ir/u, v = 3, 4, . . .) and twists (9 = 7) treated there. When comparing 
with the literature, it should be noted that our definition of the Hamiltonian by 
differentiation ([35|) rescales the Fermi- velocity to 1. 

It is also interesting to compare to results obtained by the bosonisation method. 
From e.g. [3C] one can find the excitation spectrum for the untwisted chain. The 
result for a twisted chain follows from observing that the twist boundary conditions 
on the chain modify the boundary conditions of the dual boson only. This leads to 
an energy spectrum which is exactly of the form (fl3|). 

^bosonisation = ^ + 2rf ^ + ^ E + ■ ( 4 ?) 



n=l 



Here R = is the compactification radius of the boson, m is the quantized 

momentum of the boson zero-mode, and m^ ,R are occupation numbers of the n:th 
left and right moving oscillator modes of the boson. S is the quantized dual boson 
zero-mode momentum, which coincides with our S. Comparing to fl43|) , we see that 
Bethe Ansatz strings and holes are mixtures of bosonic zero and oscillatory modes. 

4 The Kac table in Liouville theory 

As explained in the introduction, imaginary coupling Liouville theory can be 
described by a collection of twisted XXZ chains, if j/tt = \ij{y + 1) is rational 
[|17||.p| The twists are quantized in terms of 7; 

6 = kj, k — 0, ... ,u . 



2 More exactly, in \ V7 the Liouville equivalence was proved only for 7r/7 = v + 1, i^ = 2,3, 



The equivalence can be easily extended to more general rational values. 



This collection of twisted XXZ1/2 chains is exactly the one that describes the 



Bethe Ansatz of a root-of-unity Uq^sl?) invariant periodic chain [|13|, [14], [H|. This 
equivalence is not surprising, if one keeps the connection to sine-Gordon theory in 
mind. The twist-sectors in the lattice Liouville model correspond exactly to the 
#-vacua of its massive perturbation, the sine-Gordon model. These vacua become 
non-degenerate in finite volume (see e.g. [16|])- The sine-Gordon model has an 
U q {sl2) symmetry |31], fffi] , where q is related to the coupling constant [3 2 = 8(71" — 7). 

On the other hand, the integrable lattice discretation of the sine-Gordon model 
is an inhomogeneous XXZ^ chain with anisotropy 7. The massless limit, i.e. the 
imaginary coupling Liouville model, is then described by an homogeneous XXZx/2 
chain with U q (sl2) symmetry, of the type described in Refs. |TB|, [Lj|, |H) |. 



The consistent quantum group restriction of the Hilbert space of these root-of- 
unity U q (sl2) invariant models comes about by restricting k to be nonvanishing. In 
sine-Gordon theory this leads to the restricted sine-Gordon model, which flows to 
minimal models in the ultraviolet [31, Here, we are investigating the Bethe 



Ansatz description of this ultraviolet limit, and we shall see that the primary 
operators creating string states furnish the whole Kac table. 

The ground state of the quantum group restricted model is in the k = 1 sector, 
reproducing 

c = 1 - 67 2 /7r(7r - 7) . (48) 

For rational 7/7?, this gives exactly the central charges of minimal models. 

Here it should be noticed that in the quantum group invariant chain, the twist 
does not arise from boundary conditions on the spatial boundaries of the chain. 
Rather, different twists select different sectors in the Hilbert space transforming in 
a specific way under the global zZ-symmetry of the theory. In a bosonized language 
this symmetry is — > <fi + 2mr, n G 2Z. This is a symmetry of all theories with 
the potential a function of e 1 ^, including the sine-Gordon and imaginary Liouville 
ones. Thus the situation for a restricted quantum group invariant chain (i.e. a 
q-restricted Liouville chain) differs from the XXZ1/2 chain with twisted spatial 
boundary conditions discussed in the previous Section. The equation of Refs. 0, |3|], 
connecting the ground state energy to the central charge, does not get boundary 
corrections, and the central charge is as above. 

Conformal primary and secondary operators create excited states from the 
vacuum. Following Cardy [|J, the critical indices of the operators creating excited 
states are given by the finite size energies and momenta. 

We are interested in primary states with equal holomorphic and anti-holomorphic 
confomal weights, i.e. with vanishing finite size correction to the momentum 



conformal spin). Inspired by Ref. fll8 |, we look for these among the excited states 



with one higher string (n& = 1). From equation we see that a sufficient condition 
for the vanishing of the conformal spin is S = and 

kAh + I^ = . (49) 

The corresponding energies are 

1 f K 7 - 7r/ (fc)± /A; T tt(A; - l)f- -) , (50) 




2N { vr - 7 



where the sign refers to the position of the string close to A = ±00. The harmonic 
oscillator pieces in the energies vanish for all of these states. 

For generic values of 9 and 7 < tt/2, the fc-string may be close to 00 or —00, 
depending on the values of k, 7 and 9. Not all string legths are allowed for all 7 
and 9, nor all combinations of different hole and string positions. Apart from the 
Takahashi- Suzuki conditions |33[] that restrict allowed string lengths, one has to 
treat the Bethe Ansatz equation for to find the allowed values of 

For our purposes, it suffices to analyze the asymptotics of Zk(() for the string 
position ( close to ±00, which gives the possible range of quantum numbers I^ h \ 
We get different branches of primary states solving the condition ([49|). 

The first branch, for which 1^ = 0, exists when 

k + K < - . (51) 

7 

For states with hole distributions symmetric up to the twist effects of Equation (|49|) , 
this /c-string is close to —00. 

Using the formula Exited — E a = 77 (A + A) of Ref. [Q, we get the conformal 
weights of these primary states: 

x Mk - 1) - K7) 2 - 7 2 , . 

A = A = ^ '-. 4 - . 52 

47r(7r — 7) 

For unitary minimal models (7 = ^j), these states consisting of one fc-string with 
JW = and a symmetric distribution of 2k — 2 holes, furnish the whole Kac table. 
In the q-restricted case, the remnant k gets the values k — 1, . . . , v. The Takahashi- 



Suzuki condition [33] on the string length allows all strings up to k = v. For an 
unitary model, Condition (|51|) restricts this further to A; < v + 1 — k. Thus the 
integers k and k give exactly the two integer labels of the full Kac table, with the 
apropriate ranges. 

For non-unitary minimal models, these states yield only a part of the table. Some 
primaries are left out due to the technical restriction adopted in this paper, kj < ir. 

5 Conclusions 

We have computed the finite size corrections to the energies and momenta of excited 
states in XXZx/ 2 chains periodic up to a twist, for generic anisotropy 7 < 7r/2 and 
twist 9. The excited states were described as holes in the sea of real Bethe Ansatz 
rapidities, as well as complex rapidities collected in (positive parity) strings. We 
restricted outselves to the case where there are strings of only one spcified length k 
in addition to the real rapidities. 

For extracting the scaling information, we used the Destri-de Vega approach to 
treat Bethe Ansatz equations, generalizing it to cope with string-like excited states. 

Apart from completing the picture of excitation energies to be found from the 
literature, to arbitrary coupling, twist and string length (up to k < n/j), our result 



(f43l , [45"D is interesting due to its special dependence of the "field theoretic" quantum 
numbers h ± . This we consider the major novel feature of our finite size solution. 

One of the original motivations to develop the DdV formalism was the possibility 
to derive thermodynamic Bethe Ansatz (TBA) equations directly from the albegraic 



Bethe Ansatz [p7 ]. The field theoretic TBA degrees of freedom correspond to holes 
and strings, which are excitations above the vacuum of interacting magnons. The 
holes correspond to right and left moving kinks, and the strings correspond in this 
case to weakly bound kinks, with binding energy of the order ~ 1/N. 

As can be seen from Equation (f44|) , the twist dependence of the energy for a 
fixed field-theoretical configuration is not innocent. Due to the integer part cuts, 
the energies and momenta have discontinuities at some specific values of the twist. 
From the point of view of the original magnon degrees of freedom Ay of the Bethe 
Ansatz equations (|H]), however, an adiabatic excursion in 9 is possible. The A^:s 
are the building blocks of the underlying interacting vacuum, and the energies and 
momenta are continuous if the A^:s move continuously. 

During such an adiabatic excursion, the number of left and/or right holes change 
when 9 moves over a value giving a jump in the integer part expressions of (54). 
Thus at these values of 9 the vacuum absorbs a left kink state and emits a right 
one, or vice versa. It should be stressed that this spectral flow is not just a trivial 
consequence of the choice of branch in the counting functions. Rather, it encodes 
drastic changes in the structure of the Hilbert space for different values of 9. For 
example, the degeneracies of secondary states built by driving some of the holes into 
the bulk change at the jump points. 

Adiabatic excursions of this type have been investigated in the context of a XXZ 
chain threaded by a magnetic flux |2^, |24], [25]. For this, the twist is given the physical 
interpretation as a threading magnetic flux. For full comparison with the numeric 
excited state results of [|4|, |25|, negative parity strings have to be encorporated into 
our treatment, as these play a central role for larger values of 7. 

Finally, we showed how in the q-restricted Hilbert space of the U q (sl<i) invariant 
periodic chain, which is equivalent to the imaginary coupling Liouville model and 
the UV-limit of sine-Gordon, the central charge decreases from c = 1 to those of 
minimal models. Excited states with one higher string and a minimal amount of 
holes, give rise to all the primary states in the case of unitary minimal models. 

To conclude, we want to comment that the string picture should not be necessary 
for deriving the results of this paper. In Ref. p8f , wide and close pairs were used 
to describe complex excitations. This approach could certainly be generalized to 
twisted boundary conditions in the XXZ chain. The strings, however, provide a 
good tool which singles out very specific combinations of wide and close pairs that 
correspond to primary states in the critical theory. 



Acknowledgements 

We thank L.D. Faddeev, H.-P. Eckle, A. Luther, C. Montonen, A. Neveu, H.J. de 
Vega and ALB. Zamolodchikov for conversations on the subject of this paper. T.N. 
acknowledges the hospitality of NORDITA during a visit where part of this work 
was done. 



References 

[1] J. Kogut, Rev. Mod. Phys. 51 (1979) 659. 

[2] H.W. Blote, J.L. Cardy and M.P. Nightingale, Phys. Rev. Lett. 56 (1986) 742. 

[3] I. Affleck, Phys. Rev. Lett. 56 (1986) 746. 

[4] J.L. Cardy, Nucl. Phys. B 270 [FS16] (1986) 186. 

[5] A.G. Izergin and V.E. Korepin, Lett. Math. Phys. 5 (1981) 199; Nucl. Phys. 
B205 [FS5] (1982) 401. 

[6] H.J. de Vega and M. Karowski, Nucl. Phys. B 285 [FS19] (1987) 619. 

[7] F. Woynarowich and H.P Eckle, J. Phys. A 20 (1987) L97; F. Woynarowich, 
Phys. Rev. Lett. 59 (1987) 259. 

[8] N.M. Bogoliubov, A.G. Izergin and N.Yu. Reshetikhin, J. Phys. A 20 (1987) 
5361. 

[9] F.C. Alcaraz, M.N. Barber and M.T. Batchelor, Phys. Rev. Lett. 58 (1987) 
771; Ann. Phys. 182 (1988) 280. 

[10] C.J. Hamer, G.R.W. Quispel and M.T. Batchelor, J. Phys. A 20(1987) 5677. 

[11] F.C. Alcaraz, M. Baake, U. Grimm and V. Rittenberg, J. Phys. A 22 (1989) 
L5. 

[12] V. Pasquier and H. Saleur, Nucl. Phys. B 330 (1990) 523. 

[13] M. Karowski and A. Zapletal, Nucl. Phys. B 419 [FS] (1994) 567. 

[14] H. Grosse, S. Pallua, P. Prester and E. Raschhofer J. Phys. A 27 (1994) 4761. 

[15] S. Pallua and P. Prester J. Phys. A 29 (1996) 1187. 

[16] D. Fioravanti, A. Mariottini, E. Quattrini and F. Ravanini, Phys. Lett. B 390 
(1997) 243. 

[17] L.D. Faddeev and O. Tirkkonen, Nucl. Phys. B 453 (1995) 647. 

[18] M. Karowski, Nucl. Phys. B300 [FS22] (1988) 473. 

[19] H.J. de Vega and F. Woynarovich, Nucl. Phys B251 [FS13] (1985) 439. 

[20] F. Woynarovich, H.-P. Eckle and T. Truong, J. Phys. A 22 (1989) 4027. 

[21] H.-P. Eckle and C.J. Hamer, J. Phys. A 24 (1991) 191. 

[22] A. Klumper, T. Wehner and J. Zittarz, J. Phys. A 26 (1993) 2815. 

[23] B. Sutherland and B.S. Shastry, Phys. Rev. Lett. 65 (1990) 1833. 

[24] N. Yu and M. Fowler, Phys. Rev. B 46 (1992) 14583. 

[25] R.A. Rdmer, H.-P. Eckle and B. Sutherland, Phys. Rev. B 52 (1995) 1656. 

[26] P.A. Pearce and A. Klumper, Phys. Rev. Lett. 66 (1991) 974; Physica A183 
(1992) 304. 

[27] C. Destri and H.J. de Vega, Phys. Rev. Lett. 69 (1992) 2313; Nucl. Phys. 
B 438 (1995) 413. 

[28] C. Destri and H.J. de Vega, Nucl. Phys. B 504 (1997) 621. 



[29] V.E. Korepin, N.M. Bogolyubov and A.G. Izergin, Quantum Inverse Scattering 
Method and Correlation Functions (Cambridge Univ. Press, Cambridge, 1993). 

[30] S. Eggert and I. Affleck, Phys. Rev. B 46 (1992) 10866. 

[31] F. Smirnov, Nucl. Phys. B 438 (1990) 156. 

[32] D. Bernanrd and A. LeClair, Nucl. Phys. B 340 (1990) 721. 

[33] M. Takahashi and M. Suzuki, Prog. Theor. Phys. 48 (1972) 2187. 



